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METHOD FOR COMPUTER-ASSISTED DETERMINATION OF AN 
OPTIMUM-FUEL CONTROL OF NOZZLES 

CROSS-REFERENCE TO RELATED APPLICATIONS 
[0001] The present application claims priority under 35 U.S.C. §1 19 of German Patent 
Application No. 103 11 1192, filed on March 18, 2003, the disclosure of which is 
expressly incorporated by reference herein in its entirety. 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

[0002] The present invention relates to a method for the computer-assisted 
determination of an optimum-fuel control of nozzles. Such a method can be used in 
particular for controlling the nozzles of a spacecraft, such as, e.g., a satellite, a space 
probe, a space station or the like. 

2. Discussion of Background Information 

[0003] A method for the optimum-fuel control of the attitude and nutation of a spinning 
spacecraft is known from US 6,347,262. An error signal is hereby determined via a 
determination of spinning rates and spinning angles and a torque applied to the spacecraft 
according to the result for this error signal. 

[0004] EP 0 977 687 Bl describes different possibilities for the optimum-fuel 
computer-assisted control of nozzles of a spacecraft. Mainly methods are described 
thereby that contain a simplex algorithm, whereby on the other hand it is already stated 
there that such a simplex algorithm is associated with a high expenditure of computation 
time. The method described within the scope of that invention uses a method that 
depends on the simplex algorithm and starts from the formation of a simplex table. A 
dual simplex algorithm is ultimately used there to form an optimum-fuel control vector. 
As an alternative possibility for determining an optimum-fuel control of nozzles, only a 
"table look-up" method is described in EP 0 977 687 Bl, in which optimum-fuel nozzle 
arrangements are calculated and entered in a table, and a current control result is formed 
from these previously stored results for the respective current control case through the 
combination of stored results. However, it is disadvantageous that in general an actual 
optimum-fuel result is not found with this method. 
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SUMMARY OF THE INVENTION 
[0005] It is therefore the object of the present invention to provide a method for the 
computer-assisted determination of an optimum-fuel control of nozzles that manages with 
the lowest possible computation expenditure and nevertheless reliably leads to an 
optimum-fuel solution. 

[0006] This object is attained through the features of claim 1. The invention further 
comprises a computer program according to claim 5 and a computer program product 
according to claim 6. 

[0007] The present invention comprises a method for the computer-assisted 
determination of an optimum-fuel control of nozzles according to a control instruction b 
= Ax, whereby 

b represents a desired m-dimensional forces/torque vector, 
A represents an mxn-dimensional nozzle matrix and 

x represents the sought n-dimensional nozzle control vector with the aid of which 
a corresponding thrust is effected through the associated controlled nozzles. The nozzle 
control vector should thereby meet the minimization criterion 

J := E-^ -> min 

This minimization criterion is the term to indicate that the fuel consumption should be 
minimal, which is achieved through a minimization of the nozzle control and thus also a 
minimization of the thrust generated overall through the respective nozzles. 
[0008] According to the invention, the following process steps are now provided in the 
scope of this method: 

• A defined matrix transformation of starting constraints for the mass flow of the 
nozzles and of the minimization criterion takes place in a computer-assisted 
manner, whereby the starting constraints and the minimization criterion are 
respectively subjected to an analogous matrix transformation, 

• A data processing representation of a geometric description of the matrix- 
transformed starting constraints takes place in a computer-assisted manner, 
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• Through a computer-assisted geometric search procedure in the vector space, a 
computer-assisted determination of limiting point sets of the geometric 
description of the starting constraints takes place, and 

• The matrix-transformed minimization criterion is applied to the points of the 
limiting point sets. 



[0009] A simplex method or a comparable iterative method can thus be avoided within 
the scope of the invention. Instead of an iterative approach, a geometric approach 
realized in a computer-assisted manner is selected. Through the defined matrix 
transformation, a geometric description of the problem is possible that permits a 
geometric location of an optimum-fuel solution. Such a method can be much quicker 
than a solution of the problem via customary simplex methods. 
[0010] In particular it can be provided in the scope of the present invention that: 

• For the matrix transformation of the starting constraints for the mass flow of the 
nozzles, a homogenous solution of the control instruction according to Xh 0 =A 0 r is 
defined, whereby 

A Q : represents the n x (n-m) dimensional zero space matrix of A and 
r\ represents an (n-m) dimensional vector of any real numbers 

• Within the scope of the use of the matrix transformation of the minimization criterion 
a computer-assisted calculation is made of scalar products of a vector representation 
of points of the limiting point set and the vector 



v, := 



n n n 



p\-n-m 



and 

• An optimum-fuel solution is calculated with the aid of the vector r whose scalar 
product is minimal with the vector vj. 



[0011] This is an example of how through the use of a matrix transformation with the 
aid of the zero space matrix Ao, the minimization criterion as well as the starting 
constraints can be suitably transformed. The minimization criterion can thereby be 
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essentially reduced to the simple formation of scalar products of vectors, thus to a simple 
geometric calculation specification. 

[0012] A preferred further development of the present invention provides that 

• The matrix-transformed starting constraints for the mass flow of the nozzles is 
converted in a computer-assisted manner into allowable multi-dimensional value 
regions, 

• To determine the limiting point sets, a formation of at least one multi-dimensional cut 
set of the individual allowable multi-dimensional value region takes place, and 

• The limiting point sets are determined as those point sets that limit the at least one cut 
set. 

[0013] In this manner the number of the points to be considered in the scope of the 
method according to the invention can be clearly reduced, since only those points are 
considered that limit the at least one cut set. This means a further reduction of the 
computation time required for the method and thus a further advantage over the 
previously known methods. 

[0014] Advantageously the above-mentioned method can be still further simplified and 
thus the necessary computation time can be further reduced in that 

• A repeated projection of the allowable multi-dimensional value regions of the 
dimension p is made on a dimension p-1, until a projection of the allowable value 
regions on limiting intervals of a dimension p=l has been achieved and 

• Subsequently a computer-assisted search procedure carries out a computer-assisted 
determination of limiting point sets as cut set of limiting intervals. 

[0015] Through the repeated projection provided here of a dimension p on a dimension 
p-1, lower and lower dimensions (p-1, p-2, p-3, etc.) are thus gradually reached in which 
the value regions in question are represented, and thus problems that can be handled more 
easily in computational terms than they would be in the starting dimension p. In principle 
when a suitable adequately reduced dimension p>pi>l is reached, the determination of 
the limiting point sets could already take place. However, the repeated projection is 
preferably carried out until a description of the dimension p=l has been reached. Here 
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the problem can be solved in the easiest manner and with the lowest expenditure of 
computation time. 

[0016] Another subject of the current invention is a computer program for the 

computer-assisted determination of an optimum-fuel control of nozzles according to a 

control instruction b = A x , whereby: 

b represents a desired m-dimensional forces/torque vector, 

A represents an mxn-dimensional nozzle matrix and 

x represents the sought n-dimensional nozzle control vector and 

the nozzle control vector should satisfy the minimization criterion 

J := V — > min 

[0017] According to the invention it is provided that the computer program contains the 
following: 

• A first program routine for the computer-assisted execution of a defined matrix 
transformation of starting constraints for the mass flow of the nozzles and the 
minimization criterion, 

• A second program routine for the computer-assisted execution of a data processing 
representation of a geometric description of the matrix-transformed starting 
constraints, 

• A third program routine for the computer-assisted execution of a geometric search 
procedure in the vector space for the computer-assisted determination of limiting 
point sets of the geometric description of the starting constraints, 

• A fourth program routine for the computer-assisted application of the matrix- 
transformed minimization criterion to the points of the limiting point sets. 

[0018] Such a computer program is suitable for executing the above-mentioned method 
according to the invention. Other program routines can also be provided within the scope 
of this computer program which are suitable for executing one or more of the above- 
mentioned further developments of the method according to the invention. 
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A further subject of the present invention is a computer program product containing a 
machine-readable program carrier on which an above-described computer program is 
stored in the form of electronically readable control signals. The control signals can be 
stored in any suitable form, the electronic readout can then take place accordingly 
through electrical, magnetic, electromagnetic, electro-optic, or other electronic methods. 
Examples of such program carriers are magnetic tapes, diskettes, hard disks, CD-ROM or 
semiconductor components. 

[0019] Other exemplary embodiments and advantages of the present invention may be 
ascertained by reviewing the present disclosure and the accompanying drawing. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0020] The present invention is further described in the detailed description which 
follows, in reference to the noted plurality of drawings by way of non-limiting examples 
of exemplary embodiments of the present invention, in which like reference numerals 
represent similar parts throughout the several views of the drawings, and wherein: 
Fig. 1: Allowable region for the vectors r in a space of the dimension p=2, 

whereby an inequation of constraint n] (r - k i n i ) > 0 is met; 

Fig. 2: Allowable region for r in two-dimensional space taking into account the 

starting constraints for several nozzles; 

Fig. 3: Allowable region for r in the one-dimensional space taking into account 

the starting constraints for several nozzles; 

Fig. 4: Cut set of two (p-l)-dimensional planes; 

Fig. 5: Linear coordinates transformation of unit vectors; and 

Fig. 6: Comparison between computation time of a conventional method 

according to a simplex algorithm and the method according to the 
invention. 
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DETAILED DESCRIPTION OF THE PRESENT INVENTION 
[0021] The particulars shown herein are by way of example and for purposes of 
illustrative discussion of the embodiments of the present invention only and are presented 
in the cause of providing what is believed to be the most useful and readily understood 
description of the principles and conceptual aspects of the present invention. In this 
regard, no attempt is made to show structural details of the present invention in more 
detail than is necessary for the fundamental understanding of the present invention, the 
description taken with the drawings making apparent to those skilled in the art how the 
several forms of the present invention may be embodied in practice. 
[0022] The invention relates to the following problem: desired forces and/or torques are 
to be applied to a spacecraft through nozzles of the spacecraft. How should the nozzles 
be controlled so that the desired forces and/or torques can be achieved with a minimal 
amount of fuel for the nozzles? In addition the thrust of the nozzles has to lie in an 
allowable value region, thus between a minimum possible and a maximum possible 
value. Accordingly, the control values must therefore also lie in an allowable value 
region. 

[0023] This problem is a so-called linear optimization problem. Hitherto simplex 
algorithms were mostly used to solve such problems, as described at the outset, but they 
also exhibit the disadvantages described likewise at the outset. The present invention 
represents an improved possibility for solving the linear optimization problem. 
[0024] The linear optimization problem for an optimum-fuel nozzle control is based on 
the equation: 

b=Ax (2.1) 
whereby 

b represents a desired m-dimensional forces/torque vector, 

A represents an mxn-dimensional nozzle matrix and 

x represents the sought n-dimensional nozzle control vector. 

[0025] A minimization criterion is to be met for the sought unknown vector x, through 
which an optimum- fuel control is guaranteed, namely 
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J := YZ x i -> min ( 2 - 2 ) 

[0026] Furthermore, the values for the vector x must be in an allowable value region 
analogous to the allowable value region for the thrust of the nozzles: 

0<x, <\for i = l,---n. (2.3) 

[0027] Without restriction of the generality of the problem shown in (2.1), it can be 
assumed that the number of the nozzles n is greater than the number m of the desired 
forces and/or torques. 

n>m (2.4) 
and that the nozzle matrix A has full rank, 

rank(A)=m (2.5) 

[0028] An infringement of (2.5) would mean for example if all the nozzles pointed in 
exactly the same direction. However, this is always avoided by a logical arrangement of 
the nozzles on the satellite. 

[0029] Now the minimization criterion can be described: 

The general solution of the equation (2.1) x g can be written as follows: 

x g =x pa +Xho (2.6) 
whereby 

x pa special solution of (2. 1 ) 

Xho general solution of the homogenous equation of (2.1). 

[0030] A special solution x pa can be generated from (2.1) as follows: 

b = Ax = AA T (AA T )~ l Ax = AA t (aA t Y b = a[a t (aA t b)= A{x pa ) (2.7) 

[0031] From this the following special solution is obtained: 

x pa =A T {AA r Yb (2.8) 
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[0032] All the vectors Xh 0 form the general homogenous solution of the problem, for 
which vectors the following applies: 

Ax h0 = 0, (2.9) 



[0033] This homogenous solution can be written as follows: 



X ho - \ V 



(2.10) 



whereby 

A 0 : nx(n-m) zero space matrix of A 

r: (n-m) dimensional vector of any real numbers. 

[0034] With the aid of the equations (2.6, 2.8, 2.10) the minimization criterion from 
(2.2) can now be written as 



with 



Ojp 



7=1 



(2.11) 



(2.12) p:=n-m. 



(2.13) 



[0035] Equation (2.1 1) can now be interpreted as follows: for the permitted values 
of r the vector should be found whose scalar product becomes minimal with the 
vector v d . This results from the fact that the special solution as defined in equation 
(2.8) is no longer accessible to a further minimization. The present invention 
implements this search for the optimal vector r in a computer-assisted manner with 
the aid of a corresponding program routine of a computer program. 
[0036] Geometric description of the optimum-fuel solution: 
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The allowable values or the allowable region for r is obtained from the starting 
constraints for the mass flow of the nozzles according to (2.3) through a matrix 
transformation. This is formally obtained by inserting (2.6, 2.8, 2.10) in 2.3) 



-x pa <A 0 r<\-x pa 



(2.14) 



whereby 



1 : ^-dimensional vector in which all elements equal 1 

[0037] Now a geometric description of the starting constraints thus transformed can be 
made which in the scope of the present invention is realized in terms of data processing. 
The equation (2.14) can thereby be written with the following definitions in scalar form. 



/, := x 



pai 



T 



02 



On 



(2.15a-c) 



[0038] Thus the following is obtained: 



01 



r -- 



T a 0i 
\ a 0i a 0i J 



>0 



-a 



0i 



r--p-a 0l >0 

a oi a oi J 



for i=l,..n 



for i=l,..n. 



[0039] If now (2.16, 2.17) is divided by |o 0 , | 
possible - and with the following definitions 



(2.16) 



(2.17) 



due to constraint (2.5) this is always 
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._ 0/ 


for i=l, ...n 




(2.18) 




W 










0/ 


for i=n+l,.. 


.2n 


(2.19) 


W 






._ 1, 


fori=l,...n 




(2.20) 




Kl 






K 


._ u > 


for i=n+l,.. 


.2n 


(2.21) 




Kl 







the equations (2.16, 2.17) can be written as follows: 

n] (r - Xfa ) > 0 for i=l . .2n (2.22) 

[0040] Thus for a certain value of / the equation (2.22) can now be interpreted as a one- 
dimensional or multi-dimensional plane and thus shown accordingly in terms of data 
processing (see Fig. 1), whereby this plane is displaced by the vector Xn t with respect to 
the origin and is oriented perpendicular to the normal vector /?,. In the case of Fig. 1 this 
one-dimensional or multi-dimensional plane is represented as a straight line for reasons 
of simplicity, thus as a plane of the dimension 1 . This one-dimensional plane limits the 
range of the allowable values or the allowable vectors for r. If the sign of X t is negative, 
the zero vector is also contained in the allowable region for r ; however if the sign of A,- is 
positive, the zero vector is excluded from the allowable region. The normal vector n { 
points in the direction of the allowable region. The geometric description of the equation 
(2.22) is shown for p=2 in Fig. 1 for A,>0. It is noted that for each constraint i in (2.22) a 
counter-constraint i+n with n i+n ~nj exists (according to the starting constraints for the 
minimum and maximum mass flow of an individual nozzle), from which n pairs of 
constraints result that define allowable "area strips," as shown in Fig. 2. 
[0041] The common cut set of all allowable regions determines that region in the p- 
dimensional space in which r meets all 2n transformed starting constraints or all 
constraints from (2.22). In Fig. 2 this is shown for n=3 nozzles, thus for 6 constraints. If 
the cut set is formed from all allowable regions, a convex range results as cut set, i.e. that 
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for respectively two points that lie within the cut set, all points also lie on a straight line 
between these two points within the cut set. 

[0042] However, not only the region of the allowable values for r can be described 
geometrically, but also the optimum value for r, r oph which minimizes J in equation 
(2.11): as stated, J becomes minimal for that vector r that features the smallest scalar 
product or the smallest projection on the vector v</. In Fig. 2 this is given for the vector r 
that belongs to the vertex, marked with a circle, that lies on the straight line gi to which 
the normal vector nj belongs. The optimum point thus results as one of the vertices of the 
limiting lines of the cut set. The limiting lines thus form limiting point sets. 
[0043] The geometric description of this two-dimensional example can be directly 
applied analogously to higher dimensions p. 

[0044] If the problem is a three-dimensional problem, the n-transformed starting 
constraints can be described as n area regions between n plane pairs, whereby each plane 
pair comprises two planes parallel to one another and the area region lying between them 
represents the region for r that meets the relevant transformed starting constraints. The 
overall allowable region for r that meets all transformed starting constraints results as cut 
set of all n area regions and thus corresponds to a three-dimensional polygon. The 
optimum point also results here as one of the vertices of the limitation of the polygon 
which is now limited by limiting surfaces as limiting point sets. The optimum point is in 
turn that vertex whose vector r features the smallest scalar product with the vector v d . 
[0045] If it is a one-dimensional problem, the regions that are limited by the n 
transformed starting constraints can be described as one-dimensional intervals. The n 
transformed starting constraints are then described through limiting points. This is shown 
by way of example in Fig. 3. The brackets of the same size shown in Fig. 3 respectively 
form a pair of transformed starting constraints. The lower limits of the starting 
constraints are labeled with X\ 9 . . . h, the upper limits of the starting constraints with X4,. . . 
A,6 (equations 2.20, 2.21). In the example of Fig. 3 a cut set interval exists as allowable 
region for r which at the same time meets all transformed starting constraints. The 
optimum value for r results in turn from a point from the limiting point set, here as the 
right limiting point of the cut set interval. The vector r belonging to this point in turn has 
the smallest scalar product with the negative vector v d also shown in Fig. 3. 
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[0046] The same interpretations and descriptions are applicable to dimensions p>3. 
[0047] Computer-assisted determination of an optimum-fuel solution: 
[0048] A special example for determining an optimum-fuel solution will now be 
shown. A method is thereby executed as described in the scope of the invention, i.e., in 
particular a computer-assisted determination takes place of (p-l)-dimensional limiting 
point sets as geometric description of the transformed starting constraints in the form of 
planes. This special example now describes a preferred further development of the 
invention and contains in particular the following process steps: 

1. Computer-assisted determination of a first of the determined planes as limiting 
point set that probably contains the optimum-fuel solution. For further detail for 
finding this plane see below. 

2. Computer-assisted application of a linear coordinate transformation to all vectors 
r so that the normal vector of the first plane is subsequently parallel to the p-th 
unit vector of those unit vectors that span the vector space of the vectors r and so 
that the normal vector points in the same direction as this unit vector. 

3. Computer-assisted calculation of cut sets in the form of intersecting planes of the 
other determined planes with the first plane. The p-th component of each 
intersecting plane thereby has the same constant value as the p-th component of 
the first plane after the linear coordinate transformation. The clear geometric 
description of the intersecting planes can thus be simplified to a dimension p new = 
p-1. Thus only the first p-1 components are still needed to describe a vector 

4. Computer-assisted test: is p new different from 1? Then steps 1-3 are carried out 
again. 

Computer-assisted test: ifp new =l, then a computer-assisted test is carried out on 
whether an interval is present as cut set that represents an allowable region for r. 
If no such interval is present, the process is ended, a second plane is selected 
according to step land the process is carried out again. If such an interval is 
present, a computer-assisted multiplication is carried out of the limiting points of 
the interval with the vector v d now reduced to a vector of the dimension 1, thus to 
a scalar. The smaller value of the results of the multiplication is stored and, if the 
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process has already been conducted before for one or more other planes, 
compared with previously stored values. If the value stored last is smaller than a 
value stored previously, the process is started again according to step 1 with a 
new plane and correspondingly carried out. If the value stored last is not smaller, 
then the value stored previously represents the optimal solution. 

[0049] The individual computer-assisted steps will now be considered in more detail. 
[0050] Determination of the first of the determined planes according to step 1 : 
In principle the plane that contains the optimum-fuel solution cannot be clearly 
determined a priori. However, it can be clarified using the example of the case of p=2 
that certain fundamental statements are possible on the sought (p-l)-dimensional plane or 
straight line (in the case p=2 the sought plane thus has the dimension 1). For geometric 
reasons it is clear that the optimum-fuel solution always lies on one of the straight line 
sections that limit the cut set of all allowable regions, namely on those straight lines 
whose normal vector (which points in the direction of the allowable region) has the 
greatest scalar product with the vector v d and whose direction thus best corresponds to the 
direction of the vector v^. In Fig. 2 the normal vector nj best corresponds to the vector v d , 
so that the vector v d features the greatest scalar product with this normal vector. This 
scalar product can be calculated in a computer-assisted manner and the corresponding 
normal vector and thus the relevant (p-l)-dimensional plane can be determined in a 
computer-assisted manner. In this method with respect to Fig. 2 it is therefore assumed 
that the optimum-fuel solution lies on the straight line section that limits the cut set and is 
characterized by the normal vector n } . The approach is analogous in all the cases of this 
method. 

[0051] Thus in particular the following process steps are carried out: 

• Determination of the normal vector whose scalar product is maximum with the vector 

Determining whether the relevant straight line (or in general: (p-l)-dimensional 
plane) represents a limiting point set of the cut set (see below). If not, this straight 
line is abandoned and the straight line is selected whose normal vector features the 
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second largest scalar product with the vector v</. Then this straight line is tested again 
to see whether it represents a limiting point set. 

[0052] Coordinate transformation of the first plane: 

[0053] In order to simplify the computer-assisted determination of the intersecting 
planes which results as cut set of the individual (p-l)-dimensional planes, a new 
coordinate system is introduced. This is shown in Fig. 4. The coordinate system is 
rotated through a linear coordinate transformation such that a unit vector (here the unit 
vector with the highest index, i.e., index p is selected) coincides with the normal vector of 
the plane that was selected as the first plane in the step described above. Through this 
coordinate transformation an advantage results that the p-th coordinate of all intersecting 
planes now feature the same constant value. Fig. 4 shows two levels pj and p2 that in the 
case of Fig. 4 feature an intersecting plane gn of the dimension 1, thus a straight line as 
cut set. Let pi now be the first plane selected according to the method described above. 
Then through a coordinate transformation the unit vector e 3 is rotated by 180° around the 
vector e$+ni in order to obtain a transformed unit vector that now coincides with the 
vector ni (see Fig. 5). The third coordinate of gn then becomes a constant. This 
coordinate transformation can be expressed as: 



[0054] The transformation equation (2.23) can be generalized from three dimensions to 
a p-dimensional problem with 




(2.23) 



with 



d:= 




(2.24) 




(2.25) 



with 




(2.26) 
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whereby / is the index of the selected plane. 

[0055] Determination of the intersecting planes/cut sets: 

[0056] Equations for describing the cut sets can be determined in the form of plane 
equations. The plane equation for each plane i in the /^-dimensional space can be 
described as: 



[n n n i2 n i3 ....n ip _ l n ip 



< P -\ 
L r p J 



-X =0 



(2.27) 



[0057] The cut set of this plane with a selected plane j is described by the equation: 

(2.28) 



r p =Xj 



[0058] Thus this is obtained as a description of the cut set in the form of an intersecting 
plane: 



[n a n a n a —n^ 



' P -\ 



+ n ip X j -X i =0. 



(2.29) 



[0059] Equation (2.29) can be transcribed in the form of the equation (2.27) to give 



ru r -Xi = 0. 



(2.30) 



with 
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(2.31) 



r := 



'p-i 



(2.32) 



(2.33) 



[0060] Computer-assisted search procedure: 

[0061] The computer-assisted search procedure can now in principle contain the 
following steps: 

1 . Computer-assisted determination of the most promising first plane as start plane 
which most likely contains the optimum-fuel solution and storage of the vector 
that points from the point of origin to this plane. 

2. Computer-assisted execution of a coordinate transformation as described above 
and computer-assisted determination of intersecting planes as also described 
above. 

3. With the determined intersecting planes the process steps described above are 
then executed until the dimension is reduced to p=l. 

4. It is then determined whether an allowable interval is present (see above). If this 
is given, according to the above description a point is determined as description of 
an optimum-fuel solution and through a corresponding back-transformation of the 
coordinates the optimum-fuel solution is calculated in a computer-assisted 
manner. 
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5. Repetition of steps 1 through 4 until no further improvement is obtained. The 
solution then present represents the global optimum-fuel solution. 

[0062] Fig. 6 shows within the framework of a calculation example a comparison 
between the computation time of a conventional method Ml according to a simplex 
algorithm and the method M2 according to the invention, respectively for arrangements 
with 4, 5 , 6, 7 or 8 nozzles. As is already clear from Fig. 6a), in this example the 
computation time for the method M2 according to the invention for the cases with 4, 5 
and 6 nozzles is clearly shorter than the computation time of the simplex method Ml. 
Fig. 6b) illustrates the computation time Ml /computation time M2 ratio. The 
computation time of Ml is slower than the method according to the invention by a factor 
of 100 for 4 nozzles, by a factor of 16 for 5 nozzles and by a factor of 3 for 6 nozzles. 
Only for higher numbers of nozzles is the method M2 according to the invention slower 
than the simplex method in this example. 

[0063] If the optimal solution in the one-dimensional space is to be determined (i.e., if 
the difference between the number of nozzles and the number of the constraints equals 1), 
the above description of the method can be brought into a more compact form. Then 
equation (2.14) results as n equations for the upper limit and lower limit, thus as the 
allowable region of the transformed starting constraints 



with the ^-dimensional zero space vector ao and the unknown scalar r. 

[0064] If a component a 0i is equal to zero, §<=x pai <=\ must apply for the special 

solution, otherwise the above-mentioned constraint cannot be met. However, this case is 

excluded here, since this case in practice would mean a badly selected nozzle 

arrangement. 

[0065] The equations (2.23) are divided by the components dot (taking into account the 
sign of aot), and the result is 



-x pa <a Q r<l-x i 



pa 



(2.23) 



'mm r 



max 



(2.24) 
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with 



r ■ \~ s 

mm / 



pai 



«0 

l-x 



for a 0i )0 



pat 



a, 



far a 0/ <0 



0/ 



(2.25) 



r ■ \= < 

max/ 



l-x 



pai 



«0, 



pai 



for a 0i )0 



a, 



for « 0 ,<0 



0/ 



(2.26) 



[0066] The constraints in (2.24) can only be met at the same time if 



max(r m J<min(r ra J 



(2.27) 



[0067] The optimum solution for r , r opt , depends on the sign of the scalar v d in the 

minimization criterion (2.11): for v<* greater than zero, the minimization criterion is as 
small as possible if r is selected on the left edge of the interval (2.27), 



r opt = max(r min )for v d )0 



(2.28) 



and on the other hand for Vd smaller than zero, the minimization criterion is as small as 
possible if r is selected on the right edge of the interval (2.27) 



r =min(r max )for v d {0 



(2.29) 
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[0068] For v^=0 any value for r can be selected from the interval (2.27) in order to meet 
the minimization criterion. 

[0069] The values for a 0 and v d depend only on the position and thrust direction of the 
nozzles, but not on the current thrust demands. The case discrimination according to 
(2.25, 2.26, 2.28, 2.29) therefore needs to occur only once before the computer-assisted 
execution of the method to determine the optimal solution. 

[0070] It is noted that the foregoing examples have been provided merely for the 
purpose of explanation and are in no way to be construed as limiting of the present 
invention. While the present invention has been described with reference to an 
exemplary embodiment, it is understood that the words which have been used herein are 
words of description and illustration, rather than words of limitation. Changes may be 
made, within the purview of the appended claims, as presently stated and as amended, 
without departing from the scope and spirit of the present invention in its aspects. 
Although the present invention has been described herein with reference to particular 
means, materials and embodiments, the present invention is not intended to be limited to 
the particulars disclosed herein; rather, the present invention extends to all functionally 
equivalent structures, methods and uses, such as are within the scope of the appended 
claims. 
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